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Spin-crossover molecules having a low-spin ground state and a low lying excited high-spin state 
are promising components for molecular electronics. We theoretically examine one-dimensional spin- 
crossover chain molecules of the type of Fe 2+ triazole complexes. The existence of the additional 
low-spin/high-spin degree of freedom leads to rich behavior already in the ground state. We obtain 
the complete ground-state phase diagram, taking into account an elastic nearest-neighbor interac- 
tion, a ferromagnetic or antiferromagnetic exchange interaction between the magnetic ions, and an 
external magnetic field. Ground-state energies are calculated with high numerical precision using 
the density-matrix renormalization group (DMRG). Besides pure low-spin, high-spin, and alternat- 
ing low-spin/high-spin phases we obtain a number of periodic ground states with longer periods, 
which we discuss in detail. For example, for antiferromagnetic coupling there exists a dimer phase 
with a magnetic unit cell containing two high-spin ions forming a spin singlet and a single low-spin 
ion, which is stabilized by the energy gain for singlet formation. 

PACS numbers: 75.10.Jm, 75.50.Xx, 85.65.+h 



I. INTRODUCTION 

One of the most active fields of materials science to 
emerge in recent years is molecular electronics^^ which 
proposes to use individual molecules as electronic compo- 
nents. A related idea is to use a single quantum spin, say 
of an ion, to store information. Kahn and coworkers&2*£ 
have emphasized that spin- crossover compounds^^^^ 
(SCC's) are particularly promising for molecular memory 
devices. These compounds consist of complexes involv- 
ing transition-metal ions and organic ligands^aiiill The 
magnetic ions can be either in a low-spin (LS) or high- 
spin (HS) state, i.e., for the spin operator Si at site i the 
eigenvalues of Si • Si are 5 L s(5ls + l ) and ^hs^hs + l ) 
in the LS and HS state, respectively. The energy dif- 
ference between HS and LS states is due to the competi- 
tion between the crystal field splitting, which prefers dou- 
bly occupied d-orbitals and hence LS, and Hund's first 
rule, which favors the HS state, 7 i 8 i 9 i 10 In SCC's the LS 
state is the ground state and the HS state is at a mod- 
erate (thermal) excitation energy. SCC's show a char- 
acteristic crossover from the LS ground state to domi- 
nantly HS behavior at higher temperatures^ due to the 
higher degeneracy of the HS state. This crossover is typ- 
ically sharper than expected for noninteracting magnetic 
ions and is even replaced by a first-order transition in 
several compounds jLSaLifi Spin-crossover phenomena are 
also observed in organic radicals^ and certain inorganic 
transition-metal compounds^ 

Of the large number of known SCC's some naturally 
form one-dimensional chains, for example Fe 2+ with 4- 
R-l,2,4-triazole ligands. 4,13 Three ligands form bridges 
between two adjacent iron ions. Other SCC's consist of 
two-dimensional layers, for example TlSr^CoOsiiSii^ 

Besides possible applications? 3 ^ 5 - SCC's are also inter- 
esting from a statistical-physics point of view. Compared 
to conventional local-moment systems they introduce an 



additional Ising degree of freedom cr$ , which dcstinguishcs 
between the LS (cr, = +1) and HS (<7j = — 1) states. In 
the case of a diamagnetic LS state (5ls = 0) the low 
spins are essentially switched off. These SCC's are thus 
related to site-diluted spin models) 15 i 16 i 17 but m our case 
the presence or absence of a spin is a dynamical variable 
and not a type of quenched disorder. Also related are re- 
cent studies of magnetic models with mobile vacancies 1 ^ 
and of insulating phases of atoms with spin in optical 
lattices^ We show below that there is also a close rela- 
tion to finite antiferromagnetic spin chains. 

Antiferromagnetic spin chains have attracted a lot 
of interest since Haldane's (in the meantime firmly es- 
tablished) conjecture of a fundamental difference be- 
tween (isotropic) half-integer and integer quantum spin 
chains. 20 Among other things, the latter always show an 
excitation gap, while the former are critical. The valence- 
bond-solid model (AKLT model^I), in which each spin of 
length 5 is replaced by 25 fully symmetrized spin- 1/2 ob- 
jects that are then linked by singlet bonds between sites, 
was found to explain all main features of integer quan- 
tum spin chains. One peculiarity of the AKLT model is 
that at each end of open spin chains 5 of the spin-1/2 
objects find no singlet partner and form a free spin 5/2. 
For integer spins this leads in the AKLT model for even 
chain lengths to a [2(5/2) + 1] 2 = (5+l) 2 -fold degenerate 
ground state instead of the non-degenerate ground state 
found for periodic boundary conditions. This observa- 
tion carries over to antiferromagnetic Heisenberg chains. 
There, one finds a group of (5 + l) 2 low-lying states 
that become degenerate exponentially fast for long open 
chains. The lowest-lying of these states has total spin 0; 
above this state there follows a spin-1 triplet, etc. The 
maximum total spin in this group of states is given by 
5 and is concentrated at the edges. The lowest-lying ex- 
citation above them is the first true bulk excitation and 
corresponds to the lowest-lying state with M = 5+1. 
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This phenomenon has been observed experimentally, 22 
and generates a wealth of low-lying excitations if there 
are segments of spin chains of various lengths like in 
SCC's with Sls = 0. For odd chain lengths, the situ- 
ation is different as the lowest-lying states of the mag- 
netization sectors M < S are strictly degenerate both in 
the AKLT and Heisenberg models. However, the ground 
state of the magnetization S + 1 sector again contains the 
lowest-lying bulk excitation. This observation points to 
a special role of magnetization M — S, as will be seen 
throughout this paper. 

In the present paper we focus on the ground-state prop- 
erties of one-dimensional SCC's, for which we obtain es- 
sentially exact results. In particular, we treat the case 
of a diamagnetic LS state appropriate for Fe 2+ ions. We 
mostly consider antiferromagnetic coupling between the 
spins, which is probably the more common situation. 



II. THEORY 

We start from the Hamiltonian 23 

#o = -V Wj -B J2<n-hY,S?. (1) 

(ij) i i 

The sum over (ij) counts all nearest-neighbor bonds once 
and the eigenvalues of Sf are m, = — Si, — Si + 1 , . . . , Si, 
where 5, = Sls (Shs) f° r f» = 1 ( — !)• V describes 
an interaction that for V > (V < 0) favors homoge- 
neous (alternating) arrangements of LS and HS. At least 
in a subset of known SCC's this interaction is of elas- 
tic origin 2 ^ and can be of either sign^ We approximate 
this interaction by a nearest-neighbor term. 2Bq > de- 
scribes the energy difference between HS and LS and h is 
the physical magnetic field with g factor and Bohr mag- 
neton absorbed. Ho is diagonal in the basis of eigenstates 
of all Si • S l and Sf . In this basis 2 ^ 

Hq = - V Uio-j - B ^2 a i -h^^rrii. (2) 

(ij) i i 

For h = we reobtain the Ising-type model introduced 
by Wajnflasz and Pick2& for magnetic molecular com- 
pounds and by DoniachSI for lipidic chains. Here, each 
site can be in two states characterized by <j, like in the 
Ising model, but the states are degenerate with degen- 
eracies 2Sls + 1 an d 2Shs + LS The model can be rig- 
orously mapped onto an Ising model in a temperature- 
dependent effective field B^SLS^ and has been treated 
in the mean-field approximation and with Monte Carlo 
simulations A A related model with next-nearest-neighbor 
elastic interactions has recently been studied by Monte 
Carlo simulations and a number of stripe phases have 
been found^ The one-dimensional model suitable for 
triazole compounds has not been treated before. 

We are interested in a model with an additional ex- 
change interaction J ^ between the magnetic ions. The 



Hamiltonian Q is generalized to 

H = H - J^Si-Sj, (3) 

(ij) 

where J > ( J < 0) corresponds to a ferromagnetic (an- 
tiferromagnetic) coupling between the spins. J has not 
been measured in SCC's, but it has been determined in 
similar metal-organic complexes. \J\/ks is typically of 
the order of 10 to 20 K and is antiferromagnetic^SiSiiS 
as expected for a kinetic superexchange interaction^ In 
compounds based on Prussian blue, larger exchange in- 
teractions have been observed, leading to ferrimagnetic 
order at room temperature^ 

The full Hamiltonian H commutes with the operators 
Sj - Si. Thus the total spin at each site is a constant of 
motion. On the other hand, H does no longer commute 
with Sf . In the case of Fe 2+ , which is the most common 
magnetic ion in SCC's, we have Shs — 2 and Sls = 
and a significant simplification ensues, since any low 
spin partitions the chain into finite segments that do not 
interact magnetically. Thus there is a close relation to the 
physics of finite spin chains. In more than one dimension 
Sls = leads to a less trivial percolation problem — for 
long-range order to be present it is necessary for the high 
spins to percolate. 

In the following we restrict ourselves to Sls = 0- The 
Hamiltonian can be written as 

H = -Vj2(a i0 - i+1 - 1) - B J2(o-i - 1) 

i i 

-jJ2Si-S i+1 -hJ2Sf, (4) 

i i 

where we have added a constant so that the energy of 
the pure LS state vanishes. The Ising operators o~i all 
commute with the Hamiltonian H . Their eigenvalues are 
thus good quantum numbers and the Hilbert space is a 
direct product of subspaces for given {. . . , Ui, . . .}. 25 

In each sector {. . . ,<7j, a i+i, . . .} the system consists of 
chains of high spins separated by chains of low spins. The 
pure HS and LS states are obtained as the obvious limits. 
Since the LS chains do not contribute to the energy, the 
total energy in a sector can be written as a sum over 
the energies of HS chains of various lengths, including a 
contribution from their ends. These HS chains do not 
interact magnetically since Sls = 0. 

We are interested in the ground state and thus con- 
sider the lowest energy in each sector. The lowest energy 
in a sector can be written as the sum over the ground- 
state energies of non-interacting finite HS chains. Since 
H commutes with the total spin of each HS chain sep- 
arately, the z-components M of the total spins of the 
finite chains are good quantum numbers. Let us denote 
the lowest energy of a HS chain of length n with magnetic 
quantum number M by e°(M), where \M\ < nSns- We 
write 

e° n (M) = W + 2nB - hM + Ae° (M), (5) 
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where the first term comes from the extra energy of the 
change from HS to LS at the ends. The final term is the 
lowest eigenenergy of the finite HS Heisenberg chain with 
open boundary conditions and the Hamiltonian H n = 
—JJ2i=i Si ■ Si+i. 



III. RESULTS AND DISCUSSION 

A. Ferromagnetic Coupling 

For ferromagnetic coupling, J > 0, and magnetic field 
h > both the exchange interaction and the Zeeman 
term favor ferromagnetic alignment. The lowest-energy 
state thus has the maximum magnetic quantum num- 
ber M = nSas for each chain. For J > we have 
Ae°(nS H s) = -(n - 1) JS^ S and thus 

e° n {nS m ) = W + 2nB - nhS w - (n - 1) JS^. (6) 

Now let us consider a sector {. . . , cr,, a^i, . . .} for which 
the state consists of volume fractions p n of HS chains of 
length n. The HS chains have to be separated by at least 
one low spin. When counting this low spin with each HS 
chain, the volume fractions become (n + \)/np n . They 
must satisfy the constraint 



ST 71+1 ^ i 

V Pn < 1. 

^ ' 11 



(7) 



n=l 

The energy per site is 
e°(nS*Hs) 



/ 4F+J^ S 



+ 2£ - hS HS ~ JS^ S . (8) 



We have to solve the linear optimization problem of mini- 
mizing eo under the constraint Q . In the space of vectors 
(pi,P2, • • •) the region allowed by Eq. is a hyperpyra- 
mid with apex at zero and the other corners at points 
with p m = mj (m + 1) for one m and p n = for n ^ m. 
These are the points for which only chains of one single 
length are present and have the maximum volume frac- 
tion. This means that the finite HS chains are separated 
by single low spins. Since the allowed region is convex, 
the only possible solutions are its corners, except for spe- 
cial choices of parameters. Thus either p n = for all n 
(LS state) or p m = m/(m + 1) for one m and all other 
p n = 0. For the LS state we have eo = 0, whereas for the 
state with nonzero p n , 



4V- 2B { 



eo 



o 



2JS^ S 



hS 



HS 



1 



2B a — hSfts ^ 

(9) 



Examination shows that there are only three possible 
phases: (i) If W - 2B Q + 2JS^ S + hS m > and 
2B - hS HS ~ JS^s > or W - 2B + 2JS^ S + hS HS < 



and 2V + Bo — hS^is/2 > the ground state is the 
LS state, (ii) If W - 2B a + 2JS| S + hS HS > and 
2_B — hSns — JS^ S < the ground state has p n > and 
all other p m = 0, for n — > oo, which corresponds to the 
HS state, and the energy is eo = 2Bq — hSus—JS^g. Note 
that the HS state appears for any values of V and h for 
sufficiently large exchange interaction J. This is remini- 
cent of the exchange-induced Van-Vleck ferromagnetism 
m rare- earth compounds^ (iii) If W - 2B + 2 JS£ S + 
hSns < and 2V + B Q - hSns/2 < the ground state 
has pi — 1/2 and all other p m = 0. This corresponds to 
an alternating state of low and high spins. The energy is 
e = 2V + B - hS m /2. By using the LS/HS splitting 
Bq as our unit of energy, we obtain the phase diagram 
in V/Bq, h/Bo, and J /Bq shown as the J > part of 
Fig. El below. 



B. Antiferromagnetic Coupling 

In the case of antiferromagnetic coupling, J < 0, there 
is a competition between the exchange and Zeeman terms 
in Eq. QJ. Thus in principle finite HS chains of length n 
can occur with any magnetic quantum number M . The 
energy e°(M) of such a chain is given by Eq. We 
introduce volume fractions p n ,M of HS chains of length n 
with magnetic quantum number M. They must satisfy 
the constraint 

E — m<i. (io) 

n =l M=— nSus 

The energy per site is 

eo = $>n,Ml ^^ + 2B ). (11) 

The ground state for certain parameter values is deter- 
mined by the minimum of eo under the constraint 11011 . 
This is again a linear optimization problem. Except for 
accidental degeneracies, the minima occur at the corners 
of the allowed parameter region. Thus either p n ,M = 
for all (n, M) (LS state) or p n ,M = n/(n + 1) for one 
(n, M) and p n ,M = for all others. In the latter case the 
energy per site is 

AV-2B -hM + Ae° n (M) 

e = — h ZB . [LZ) 

n + 1 

Ae°(M)/\J\ is the lowest energy of the finite antiferro- 
magnetic Heisenberg chain with open boundary condi- 
tions and the Hamiltonian H' n = ~z~2™ = i Si ■ Sj+i in the 
sector with total S z quantum number M. It is not pos- 
sible to find these energies in analytical form. For suffi- 
ciently small n, the Hamiltonian H' n can be diagonalized 
numerically. We have calculated the energies up to n = 8 
for all M using the Lanczos algorithm. The results for 
Ae° (M)/| J\/ (n+1) are shown in Fig.Q]as colored circles, 
where identical colors denote the same values of M. 
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FIG. 1: (Color online). The lowest energies Ae°(M)/|J| 
of antiferromagnetic Heisenberg chains of length n for spin 
Shs — 2 in sectors with open boundary conditions for fixed 
total S z quantum number M. The energies are normalized 
by a factor l/(n + 1) and shown as a function of l/(n + 1), 
where n + 1 is the period of the LS/HS pattern. Circles: Re- 
sults from exact Lanczos diagonalization. Squares: Results 
from DMRG. Equal colors correspond to equal M. For odd n 
the energies for M — 0, 1, and 2 are degenerate, as noted in 
Sec. [I] The cross at l/(n + 1) =0 denotes the extrapolated 
energy density —4.761248(1) of an infinite chaim— 



The energies for longer chains can be calculated 
with excellent precision with a finite-chain DMRG al- 
gorithm(2L2& for a detailed e xplanat ion of the algorithm 
and its applications see Refs. I39l40i To obtain typically 
seven-digit precision for the ground state energy per site, 
we have kept up to M = 300 states in the reduced DMRG 
Hilbert spaces and carried out three finite-system sweeps 
which was enough to ensure convergence. Note that 
DMRG prefers open to periodic boundary conditions. In 
standard DMRG applications to integer-spin chains, end 
spins of length 5/2 (a spin 1 at each end for our case) 
are attached to eliminate the peculiar boundary degrees 
of freedom and access bulk physics directlyi^L^i In the 
present calculation, these boundary degrees of freedom 
are physical and hence no end spins are attached. The 
energies for lengths n = 9 through n — 49 as well as 
n = 99 are shown in Fig. ^ by squares. 42 

We first consider the case of vanishing magnetic 
field, h = 0. Figure ^ shows that the state with 
M = has the lowest energy for any n. Then 
e a = (4V-2B + Ae°(0))/(n + l) + 2B is a sum of 
Ae°(0)/(n + 1) and a linear function in l/(n + 1). The 
minimum of eo can only occur for n = 1, n = 2, 
or n — > oo, since all other points lie above the dot- 
ted straight lines connecting the corresponding points in 
Fig. (not obvious on this scale) . The relevant energies 
per site are thus determined by Ae^/2 = 0, Ae°/3 = 2J, 
and lim„_>oo Ae°/(n + 1) = 4.761248 J and have to be 
compared to the LS energy eo = 0. The resulting phase 
diagram is shown in Fig. [21 Note the appearance of a 



dimer (n — 2) phase. In this phase the energy increase 
due to the HS-HS pairs (V < favors HS-LS neighbors) 
is overcompensated by the large negative singlet forma- 
tion energy of Heisenberg spin pairs. 
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FIG. 2: (Color online). Ground-state phase diagram of the ID 
spin-crossover model with Sls = and Shs = 2 for vanishing 
magnetic field, h = 0, and antiferromagnetic exchange inter- 
action, J < 0. The dimer (n = 2) phase case is highlighted. 
The heavy solid lines denote discontinuous transitions. The 
various spin structures are indicated by cartoons (solid sym- 
bols: HS state, open symbols: LS state); these should not be 
overinterpreted — there is no magnetic long-range order. 

For general magnetic field h we have to take all pos- 
sible magnetic quantum numbers M of the chains into 
account. This is obviously impossible for the pure HS 
phase. Instead, we have performed DMRG calculations 
for chain length n = 99 for all possible magnetizations 
M = 0, . . . , 198 and use them as a caricature of the HS 
state. The resulting errors are discussed below. 

For each set of parameters (V/Bq, J/ Bo, h/Bg) we cal- 
culate the minimum energy densities for all states with 
n < 49 as well as n — 99 from Eq. I|12l) . The energy 
density of the LS state is zero. Then the ground state 
is obtained by finding the minimum energy. Figure 
shows a series of phase diagrams for fixed exchange in- 
teraction J. Note that the lower edges of each diagram, 
i.e., h = 0, are consistent with Fig. We observe that 
the dimer (n = 2) phase present at h = is suppressed 
by the field, as is expected since this phase is stabilized 
by the singlet formation energy. 

For J < —0.6 the phase diagrams remain qualitatively 
the same. The features are shifted to lower V and ex- 
panded linearly in both the V and h directions. Letting 
V, J, and h go to infinity while keeping their ratios fixed 
corresponds to the limit B$ — > 0, i.e., vanishing energy 
difference between LS and HS states. In this limit we 
choose | J | as our unit of energy, leaving two dimension- 
less parameters V/\J\ and h/\J\. The resulting phase 
diagram is shown in Fig. 

Interesting behavior is seen in the triangular region 
surrounded by phases with n — 1, n — 2, and the HS 
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FIG. 3: Zero-temperature phase diagrams for the same model 
as in Fig. but in a magnetic field, for antiferromagnetic ex- 
change interactions J/B = 0.0, -0.2, -0.4, -0.41, -0.42, 
—0.6. The white area corresponds to the LS phase, the black 
to the HS phase, approximated by a phase with n — 99, and 
the gray areas correspond to n = 1, 2, 3, 5 (from light to dark). 
All transition are discontinuous, the purely magnetic contin- 
uous transition discussed below is not shown. 
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FIG. 4: (a) Zero-temperature phase diagrams as in Fig. [3 in 
the limit Bo — > (or V,J,h — > oo with their ratios fixed). 
The gray scale is the same as in Fig. [3] (b) Enlargement of 
the left figure on a different gray scale. The values of n in the 
various ground states are indicated. 



phase, as shown in Fig.^Jb). Here, phases with HS chain 
lengths n = 3, 5, 7, and 9 are found. We do not observe 
any further phases. To understand why only odd n ap- 
pear, we refer to the energy densities in Fig.^ The en- 
ergy of the singlet (M = 0) state is lower for even n than 
expected from a linear fit, at least for the small n relevant 
here, while for odd n it is higher. Thus at zero magnetic 
field, even-n states are preferred. On the other hand, the 
energy of states with M > 2 is higher for even n than a 
linear fit, while for odd n it is lower. Thus in a sufficiently 
large magnetic field odd-n states are preferred. The fact 
that the series of odd n is cut off at n = 9 is a result of 
the detailed numerical values of energies in Fig. for 
larger n, the HS state happens to have the lower energy. 
While we expect the appearance of only odd HS chain 
lengths to be a robust feature of spin-crossover chains 
with strong antiferromagnetic interaction, the restriction 



to n < 9 should thus be model dependent. For example, 
inclusion of a next-nearest-neighbor elastic interaction or 
a different integer value of Shs may change this result. 
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FIG. 5: Left: Density plot of the magnetic quantum number 
M of the finite chains in the same parameter region as in 
Ffg-EI a )- Black corresponds to M — 0, white to maximum M. 
The phases with odd n all have M — 2. Right: Magnetization 
m = M/(n + 1) derived from the data in the left plot. Black 
(white) corresponds to m = (m = 1). The magnetization 
of the fully polarized HS state, which does not appear in the 
plot, would be m = Shs = 2. 



C. Magnetic Properties 

We now discuss the magnetic properties in more detail. 
In the ferromagnetic case all high spins are fully aligned 
with a nonzero magnetic field h. For the antiferromag- 
netic case Fig.[S]shows the magnetic quantum number M 
of the finite HS chains and the magnetization m in the 
ground states in the limit of large V, J, h. The magne- 
tization is defined as the magnetic quantum number M 
divided by the period, m = M/(n + 1). Interestingly, the 
phases with odd n all have M = Shs = 2, including the 
n = 7 and n = 9 phases not resolved in Fig. [3] This of 
course corresponds to different magnetizations m. To un- 
derstand the special significance of the value M = Shs, 
we plot in Fig.^lthe local expectation values of the spins, 
(Sf), for each site of a HS chain of length n = 9, obtained 
with DMRG. The plot shows that for < M < S H s the 
odd chain can accomodate the finite spin by forming a 
Neel-type state. For higher M this is no longer possible 
and spins pointing in the "wrong" direction are reduced 
(a bulk magnon is excited). Due to the cost in exchange 
energy such states are always higher in energy than com- 
peting phases. Compare also the discussion in Sec. [I] 

The dimer (n = 2) phase always consists of singlets, 
AI = 0. This shows that it is energetically favorable to 
replace the dimer state by a state with odd n or the HS 
state, instead of having M > for the dimers. 

Finally, we turn to the pure HS phase. At T = 
the system is equivalent to an infinite antiferromagnetic 
S = Shs chain. We thus expect the magnetization to 
rise continuously with increasing magnetic field h and to 
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FIG. 6: (Color online). Local expectation values (Sf) from 
DMRG for each site of a HS chain of length n — 9 for various 
total magnetic quantum numbers of the chain, M. 



reach its maximum value m = ShSi i- e - ; full spin align- 
ment, at a continuous phase transition. Since we approx- 
imate the HS phase by the n — 99 phase, the continuous 
increase is replaced by small steps. The position of this 
transition is determined by equating the energies per site 
for M = nSas and M = nS m - 1- From Eq. (H2J we 
thus obtain the critical field h c = Ae° (2n) — Ae° (2n — 1), 
which is proportional to J and independent of Bq and V. 
For n = 99 exact diagonalization yields h c « —7.9980 J, 
compared to the exact result for an infinite chain, h c = 
— 4 JShs = —8 J. To find the critical behavior close to 
this transition we define the deviation of the magnetiza- 
tion from its maximum by Am = Sus—m. Plotting Am 2 
vs. h (not shown) we find that Am? is linear in h c — h 
so that the critical exponent of Am with respect to the 
field h is mean-field-like, = 1/2. 

The previous discussion shows that by restricting the 
DMRG calculations to n < 100 we make an error for the 
transition to full spin alignment of the order of 0.03%. As 
another way to estimate the errors, we have determined 
the triple point between LS, HS, and dimer phases in zero 
field and compared the result to the "exact" triple point 
shown in Fig. [21 The error is of the order of 0.2% for V 
and 0.1% for J. 

We also obtain spin correlations from the DMRG. Fig- 
urcQshows spin-spin correlation functions for n = 99 for 
two spins close to the center of the chain, where the infi- 
nite chain should be well approximated. We first notice 
the anomaly at M = 2. This is of the same origin as the 
stabilization of M = 2 for small odd chain lengths, dis- 
cussed above. It is thus a finite-size effect not present for 
the true HS phase. Apart from this anomaly, the trans- 
verse correlations (Sf Sj~) first grow with magnetic field 
h or magnetization. This is the one-dimensional analog 
of the spin-flop state in ordered antiferromagnets, where 
the staggered magnetization is oriented perpendicularly 
to the applied field. For large fields, the correlations de- 
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FIG. 7: (Color online). Expectation value of the z-component 
of the spins, (S z ), and spin-spin correlation functions at the 
separation of \i — j\ = 4 as functions of the total magnetic 
quantum number M for chain length n = 99. Inset: Corre- 
lation function (Sf S~) as a function of separation for three 
values of the total magnetic quantum number M. 

crease again since the spins are more and more forced 
into the field direction. In the fully polarized state for 
h > h c the transverse fluctuations vanish. 



IV. SUMMARY AND CONCLUSIONS 

To conclude, we have studied a model for spin-cross- 
over compounds forming one-dimensional chains. We 
consider spin quantum numbers appropriate for Fe 2+ 
ions, which have a spin-0 LS and a spin-2 HS state. The 
model includes elastic and exchange interactions and an 
applied magnetic field. The most important effect left out 
here is probably the dipole-dipole interaction, which is of 
long range for an isolated chain, but becomes screened 
if the chain is deposited on a conducting substrate. We 
obtain the ground-state phase diagram analytically for 
ferromagnetic or zero exchange interaction and using the 
density-matrix renormalization group (DMRG) for anti- 
ferromagnetic exchange. As a summary, Fig. [S] shows 
the full phase diagram. The continuous transition to 
full spin alignment is indicated by the solid blue surface. 
All other surfaces are discontinuous transitions between 
states with different chain length n. Horizontal cuts cor- 
respond to the plots in Fig. |3 the vertical cut at h = 
to Fig. |3 Besides a diamagnetic LS phase and a HS 
phase equivalent to the usual Heisenberg chain we find 
a number of more complex phases. For sufficiently neg- 
ative elastic interaction V we find an alternating phase 
of low and high spins. In quasi- two-dimensional SCC's 
the corresponding checkerboard state has been observed 
experimentally. 12 

For antiferromagnetic coupling we find a robust dimer 
(n = 2) phase, which consists of spin singlets formed by 
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FIG. 8: (Color online). Zero-temperature phase diagram of 
the spin-crossover chain with Sls = and 5hs = 2. Posi- 
tive (negative) J corresponds to ferromagnetic (antiferromag- 
netic) exchange interaction. V denotes the elastic interaction, 
h the applied magnetic field in units of energy, and 2Bo is the 
energy difference between HS and LS states in the absence 
of interactions. The solid surfaces denote phase transitions 
between different phases, which are indicated. The phases 
with HS chain lengths n = 3, 5, . . . are hidden in this view. 
All transitions are discontinuous, except for the continuous 
transition to full spin alignment in the HS phase, shown as 
the monochrome (blue) surface. 



two high spins separated by single low spins. Since this 
phase appears at zero and low magnetic fields, it should 
be accessible experimentally. As the magnetic field is in- 
creased, the dimers remain in the singlet state until states 
with an odd number of HS ions or the pure HS state 
become lower in energy, whereupon the dimer phase is 
destroyed in a discontinuous transition. At higher mag- 
netic fields we find a number of phases consisting of fi- 
nite chains of length n — 3, 5, 7, 9 of HS ions with total 
S z quantum number M = 2 separated by single LS ions. 
We suggest that the succession of odd chain lengths is a 
general feature of spin-crossover chains. 

We thus find that a model that contains the most 
important ingredients of one-dimensional spin-crossover 
systems shows a rich ground-state phase diagram. The 
model is related to various systems studied in recent 
years, such as site-diluted spin models and finite anti- 
ferromagnetic Heisenberg chains. Questions for the fu- 
ture concern the behavior at nonzero temperature and 
of higher-dimensional models, in which percolation plays 
an important role. New physics comes into play since 
the dilution by LS ions is not quenched disorder but a 
dynamical degree of freedom. 
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